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ABSTRACT 


The thesis develops a theoretical model for predicting 
ship interactions that should be useful in a ship trajec- 
tory analysis. The method assumes that when ships are 
close together, the potential field is the main source of 
the interaction forces, The flows around the ships are 
represented by using half bodies of revolution generated 
with axial singularity distributions. A second axial 
distribution of doublets were included to correct for 
cross flows and thus maintain the rigid body boundary 
condition, following the determination of the fluid flows, 
the forces were calculated using Lagally’s theorem as ex- 
panded by W, E, Cummins, 


The principal result was a computerized model to pre- 
dict forces resulting from steady, parallel motion of 
bodies of revolution moving in an infinite ideal fluid, 
The program 1S capable of handling ship geometries in rec- 
tilinear motion. It was found that the theoretical model 
produced similar results when compared with model tests 
conducted by D, W, Taylor, however, those tests were re- 
ported as being crude, This emphasizes the need to obtain 
a set of more valid empirical data, 
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INTRODUCTION 


it General 


Interactions between ships operating in proximity can 
be thought of as a modern day phenomena, There are no knovm 
references to it in the literature until after the emer- 
gence of steam driven ships. Along with the increase in 
shipping in shallow and restricted waters, this innovation 
contributed to the higher frequency of collisions. These 
collisions provided the operator with his initial encounters 
with ship interaction forces, It was many years, however, 
before the courts recognized the existence of this strange 
force which contributed to many of the collisions. The 
acceptance of these forces came only after many eyewitness 
accounts described strange ship behaviors and it was con- 
cluded that there must have been a force that was not 
present during normal operations. 

With the increased density of shipping in restricted 
waters, especially for channels, it is not surprising that 
collisions involving interactions were first noticed for 
the shallow water occurrences. The effect of shallow water 
has persisted as being important in describing the inter- 
action phenomena. Shallow water does influence the flow 
around the ship, however, for depths greater than twice the 


draft of the ship, the effect of the ocean bottom is 
8 





expected to be small. This is based upon the interaction 
between ships in deep water forces have been experimentally 
determined as the ship separation was increased, 

With the increased speeds and maneuverability intro- 
duced by steam propulsion, interactions other than with 
ships were first noticed and recorded. These involved the 
forces experienced by the position of the ship in the 
channel. From these forces, rules of thumb were formed such 
as the ship "Smells the bank". This particular rule indi- 
cates the mariner realized that the ship interacts with the 
channel bank so as to sheer away from it. It is in this 
area of ship maneuverability in restricted water that a 
great many empirical studies have been made. With the re- 
sults of Bernoulli, the mechanism could be explained. Even 
today the operators aware of the Venturi effect tend to 
associated all interactions and explain them consistent 
with this law, This may give a rough guide, but does not 
contribute to the understanding of the process. It may at 
times even lead an investigator to the wrong conclusions. 

The forces evolving from ship interactions has long 
been a controversial subject. The main interest has stemmed 
from the cases which result from collisions. The literature 
is sprinkled with accounts of major collisions. Included 
are discussions where important figures of the day give 


conflicting explanations of the ships' behavior. The first 





scientific investigations involving model tests occurred 

in 1909 when D. W. Taylor presented his findings to the 
society of Naval Architects and Marine Engineers. The cis- 
cussions of the paper demonstrated the debate at the time 
about the causes of the forces, 

The study of interaction forces is important for at 
least four reasons, First, they contribute to collisions, 
The understanding of these forces is necessary to help pre- 
vent future occurrences. Second is for the control of 
Ships moving in proximity to each other, as done by several 
Navys to accomplish refueling of ships at sea. The third 
implication is in the design of future ships and their 
operating procedures, The adverse effects of ship inter- 
actions may be reduced by modifying design parameters and 
effective design and use of controls. Finally, for those 
ships in operation, the present rules for navigation do 
not recognize the interaction effects, Future study may 
show the need to modify these rules so as to provide safer 
movement of the merchant fleet. 

With this perspective, the research effort here at- 
tempts to determine the extent of the past efforts to de- 
velop a generalized theory for ship interactions which 
would be practical to the engineer. The view is toward a 
general model, in particular one which could be used as 
part of a generalized program for predicting ship trajecto- 


ries, 
10 





1.2 Ship Interaction 


me. . 6GeNneral 

What is meant by interaction forces? Any force that is 
only observed when one or more ships (or other objects such 
as wall, etc.) are present can be classified as an inter- 
action force. The motion of a ship through the water mani-~ 
fests itself by waves and the wake formed and by a moving 
pressure disturbance in the water. The moving pressure 
field will be identified as the potential disturbance where 
the waves will be referred to as the free surface disturb- 
ance, 

The interactions between ships operating on the surface 
of the ocean result from the influence of both the free sur- 
face and potential disturbances, whereas ‘each appears to 
create forces by different mechanisms and differ in their 
spatial distribution. The potential field generates a 
distribution of pressures in the fluid, The waves also 

generate a pressure field, However, this rapidly attenuates 
with depth. The waves also change the water surface and 
this introduces forces. If the frequency is low, the long 
duration of the force will cause significant ship motions. 
An example is in the case of following seas where direc- 
tional instability may occur and result in a large response 
of the ship. 


Very close to the ship, the potential field dominates 


al 





the fluid flow pattern, where much further away the waves 
produce the flow disturbances, The induced flows act on 
other bodies (ships) placed in the perturbed stream to 


Create interaction forces, 


to 2eie layvior Results 


One of the first scientific experiments performed and 
reported on the interaction between ships was conducted by 
David W, Taylor at the Model Basin in Washington, D, C, 
The results are reported in the 1909 Transaction of 
society of Naval Architects and Marine Engineers, A good 
summary of the work can be found in Reference 2, Here 
attention is focused on the quantitative measurements, 

The results are presented for the two model trials. The 
first 1S with relatively thin models and the second is with 
fuller models, For convenience, the model characteristics 
are also presented, In each case, the forces are reported 


for only one of the models as indicated, 
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The original results were non-dimensionalized by the 
resistance of the model, Taylor comments that the results 
may have inaccuracies due to problems encountered in towing 
the models on a parallel course, He also stated that the 
results for the thinner models were less consistant, It is 
difficult to say if the data was erratic or if the results 
were counter intuitive, thus indicating to Taylor that they 
might be erratic. Whatever the meaning, both results are 
presented and evaluated using the proposed theoretical 
models, Appendix E contains calculations of the possible 
errors in the data that might have resulted from misalign- 
ment in the models, 

The original data was modified to develop the total 
forces and moments on the models. Also a change in the 
non-Gimensionalization was made, To accomplish these 
changes, the force acting on the model was determined by 
multiplying the data by the resistance of the model. Taylor 
states that the models were towed at speed between two and 
three knots. <A value of 2,5 knots was assumed, The 
total force on the model is determined by adding the forces 
that were measured at each end of the model, The moments 
were calculated about the midship position, The following 
curves depict Taylor's data after the indicated changes, 

The curves are drawn for the ships in various rela- 


tive positions, The abscissa is scaled in terms of the 
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model length. The attraction and repulsion forces are self- 
explanatory and are applicable for the model passing on 

the starboard and port side. The moments are depicted 

as being either bow in or bow out. This indicates the di- 
rection in which the moment is acting on the model, The 
Thin Ship data represents Model Number 834 as it passes 
Model Number 838, The Full Ship data represents Model 


Number 858 passing Model Number 866, 


1.2.3 Newton Results 


Experiments were carried out at the Admiralty Experi- 
mental Works (AEW), Haslar, England, between 1946 and 1948 
to test out the feasibility of the replenishment of war- 
Ships with fuel tankers, These results were reported at 
the first Symposium on Ship Maneuverability at the David 
Taylor Model Basin in Washington, D. C. in 1960 by R. N, 
Newton. The tests conducted included an experiment with 
constrained models of the Battleship H.M.S. KING GEORGE 
and the R.F.A. OLNA. (Ship A and B Respectively) 

A.E.W. used freely propelled and controlable models to 
study the behavior during maneuvering alongside and breaking 
away. The work also contains some data resulting from 
tests conducted at sea with ships, The results are presen- 
ted here for the constrained models, The force moment act- 


ing on each ship has been converted to a position amidship. 
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FIGURE 3 


AEW MODEL TESTS 
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1.3 Previous Analytic Models 


ioe! Bl lipso1d 


Sir H. T. Havelock proposed a model of Ship interac- 
tions based upon potential flow. As an example, he cited 
the ellipsoid. In the discussion to Reference 8, the solu- 
tions were given for two unequally sized ellipsoids moving 
parallel in the fluid. The forces were calculated using 
Lagally's theorem, 

The following singularity distribution describes an 


ellipsoid in a uniform flow. 


SMe) = 2 5 Cue) 4x 
WHERE 


1/o = 4 e/(1-e“) - 2 log [ (+ e/ (1-e) } 


Vv 


Soced aay FT/SEC 


x 


Distance from midbody position 


e= Ellipsoid eccentricity 


Appendix D contains sample distributions for the 
ellipsoid model, Also, the results for the sample case 


cited by Havelock are recalculated in Appendix IF. 


18 





1.3.2 Rankimne® Ovoid 


The same model proposed by Havelock for use with the 
ellipsoid could be applied to any other body shape for 
which the singularity distributions are known, The Rankine 
Ovoid is one of the simplest bodies to generate from sin- 
gularities in the fluid. It is formed by placing two equal 
but opposite singularities in a uniform flow so that their 
axis is in line with the incoming flow and the source is 
upstream of the sink, Given the singularities, it is easy 
to determine the size of the body that is represented in 
the fluid. Silverstein in Reference 14 developed the force 
on a Rankine Ovoid. Appendix F contains the recomputation 
of a case demonstrated in his paper. Even though the dis- 
tribution is relatively simple and the body shape can be 
determined, it takes some effort to solve the inverse prob- 
lem, The desSired singularity strengths and positions for 
a given shape requires the simultaneous solution of two 
non-linear equations, 


The following equations describe the Rankine Ovoid: 


a 4, (Strength) (Dist) (L/2) 
[ >.> 
maya ~ (Dist) |@ 


es (Strength) (Dist) 
v J(B/2)2 + (Dist) 


WHERE 
V =sSpecdminett/sec- 
Dist = Distance from the modbody that singulari- 
ties are located, 


Strength = Singularity seecuet! ic 





The solution of the Simultaneous equations was accom- 
plished on the computer .using s system provided program 
called Zeroin. A copy of the instruction sheet for this 
program 1S included in Appendix H. A tabulation of Ran- 
kine Ovoid distributions 1S provided in Appendix D, 

The experience with the system routine that solved the 
Simultaneous equations should help future investigations. 
The Zeroin Program is one of two that were available. MThis 
particular one was selected mainly because it offered an 
automatic feature which permitted retrieval of intermediate 
steps in the solution so that the program results could be 
verified, It was found that an approximate solution to 
the problem using the semi-infinite body would provide a 
good starting solution for the program. The routine would 
not converge to the problem solution, however, until the 
magnitude of the initial solution vector was reduced in 
Size by two orders of magnitude. This prevented the prog- 
ram from jumping past the solution and attempting to find 


another solution at very large values of the variables. 


1.3.3 Slender Body Approximation 


Newman has derived the forces and moments acting ona 
Slender body of revolution moving parallel to a wall. In 
his derivation, it was assumed that the fluid was ideal and 


of infinite extent, thus ignoring any free surface affects. 
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In order that the boundary condition at the wall be satis- 
file, an image singularity distribution was provided. This 
is analogous to two bodies of revolution moving parallel 
and abreast through an infinite ideal fluid, 

The approach used was to determine the singularity 
distribution that would develop fluid boundaries representing 
the rigid body moving through the medium. The singularities 
were sized uSing a two dimensional model and then trans- 
formed to three dimensional by assuming that the body was 
slender (ro<< L). The location of the singularities was 
determined by satisfying the linearized boundary condition 
on the body surface. It was pointed out that for the body 
close to the wall, the singularity distribution could be 
approximated by a general bipolar distribution. This re- 
sulted in the Singularities being placed along the length 
of the body in such a way aS to be displaced from the center 
axis. This formulation requires the wall be close 
(a << 1), The displacement of the singularities is des- 
cribed by the following formulas, where the displacement 
from the centerline axis is equal to (Zo - a). 


ae = Ae e ro* 


Once the singularity distributions are known, the 
forces acting on the body are calculated using Lagally's 
theory. It 1s easy to verify that the singularity distri- 
bution remains within the body and thus is amenable to a 
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Lagally force calculation. The following diagrams show the 
distributions for a sample ship form, 

The procedure developed the singularity distribution 
using a slender body assumption, The resulting force in- 
tegrals were then modified, An approximation to the inte- 
gral is based upon the assumption that the part of the in- 
tegrand which represents the singularity as a function of 
lenzth can be treated as a constant value. 

Newman's work produced the following results for the 


problem. 


Fa 


tl 


+°/2 7! 
a TT P U e f [ y. (x) ree) | ; ; ce [row )d { : dx 
=", 


‘h 2 [ 
! ont Ao Tee 
M a. P ie in Prox) tox) J } 25 lr &)J ¢ xX dx 


yf 


y 


WHERE Lie 
I 


separation of the body axis from the wall 


Radial dimension at position X along the 
body 


Newman also presented a sample case to demonstrate the 
results, The body was assumed to be an ellipsoid des- 
cribed in the following form: 


ro (x) = ron * (1-4 rg Z) 


WHERE 


Maximum radius of the body 
Longitudinal position measured from the 


4 
tt 


midship station, 
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Substituting into the Force Equation, 


P, = - t/a TP VU ver Lb LQ. (8) - a. (8)] 

WHERE 

Q, is the Legender function of the Second Kind. 
AND &= *%, 

eae! 

ee elen - 


Fa max 
The force integrals were evaluated for this case and 


the results presented as follows, 
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In the concluding comments of the paper, the author 
states, ",...it seems likely that for bodies with full or 
blunt ends this [moment] integral will be dominated by the 
large values of S’(x) at the two ends, and thus the pitch 
moment will be positive (bow down) if the bow is blunt 
compared to the stern, and vice versa", In this comment, 
the pitch is analogous to the yaw in the previous defined 
orientations. He is predicting a 'bow out’ condition if 
the sterns are blunt compared to the bows, The author 
also notes that the force is always attractive and ".,. for 
geometrically affine bodies, having the same cross sec- 
tional distributions but different lengths, the force will 
be inversely proportional to the length whereas the moment 
will be independent of the length", "rom this result, 
the nondimensional value of maximum force is inversely 


proportional to the cube of the L/B ratio. 
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FIGURE 6 
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CHAPTER Ii. 
FORCE MODEL 


Zeit WMWodel Descripsion 


The proposed model for predicting ship interaction 
forces uses the hypothesis proposed by Havelock that the 
potential field could explain the ship interaction pheno- 
mena. This approach seems reasonable for several reasons. 
' First, the ship moving at the water's surface produces a 
flow field similar to that produced in an infinite fluid, 
Second, the potential field dominates the flow near the 
body. Pressures in this region fluctuate spatially and the 
changes find relief on the surface. The movement of this 
pressure field causes the wave system. Third, for most 
cases of interest, the ships are operating below their de- 
sign speed. This is particularly true for naval ships re- 
plenishing at sea. The fine entry of these ships causes 
minimal wave disturbances when at replenishment speeds. 

The potential flows without a free surface can be 
modeled using bodies moving in an infinite fluid with image 
systems to maintain the flow conditions on the waterline 
plane. This condition requires that there be no normai 
component of flow at the interface representing the undis- 
turbed water level, <A body of revolution moving horizon- 
tally in the fluid has a symmetry of flow and is used to 


model the ships. The maximum diameter of this body is taken 
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as the ship's beam because of the major influence this di- 
mension has. 

The force acting on the total body of revolution would 
represent twice the force experienced by the ship. The 
force results from the integration of pressures on (the 
lower) half of the body of revolution, The symmetry of 
flow will insure the same pressures and thus forces on the 
other half body. 

The major reason for selecting a body of revolution was 
uneeability of using an axial distribution of singularities 
to describe the flow, Hess et, al. have shown the compu- 
tations required for representing an arbitrary bocy in an 
ideal fluid by a surface distribution of singularities. 

The computational effort for each ship would be much greater 
for surface distribution and as such is not desirable for 
use in a Ship trajectory model, 

In total, four models were developed and programmed, 
The first two incorporated the ellipsoid and Rankine Ovoid 
as discussed earlier, The third and final one involved a 
modified ellipsoid and slender body respectively, each with 
cross flow corrections. The first two used the derivations 
of Havelock and Silverstien. The third is only a special 
case of the fourth model. The slender body model will be 
discussed in more detail, Appendix I contains a logic flow 


chart of the computer program for this model, 
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The slender body model uses bodies of revolution to 
model the ships as described previously, the axial distri- 
bution of sources and Sinks is obtained by assuming the 
body of revolution is slender. In two dimensions, the sin- 
gularities are sized using the following equations. The 


derivation for three dimensions are presented in Appendix C. 






<< UNIEORP™ 


AXIAL $ FLOW 
eee ee ARTY Se + —_____ 
——e ee See ae Se = 
DisTRIBUTION (4x) Oey X = (-V) 
6 
ee Pao 
SURFACE 
Atr ye Calais 
aia . Ox 
Xo 
WHERE 
Sve) = SOURCE STREAGTH PER UNIT LENGTH 


oe Flow MEASURED IN FT./ sec. 


This axial distribution of singularities determines the 
body in a uniform flow, When another body comes near, the 
flow field is changed, so that it is no longer uniform, 

As the bodies approach each other, their affect in altering 
the flow field causes a deformation in the fluid boundary 
which describes the rigid bodies. To handle this problem 
an additional distribution of dipole singularities are 
placed along the axis. Thes’ aire oriented in the plus or 
minus direction on the body axis and are sized to counter 
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the effects of local cross flow on the bodies, The change 
in the flow in the axial direction is handled py perturbing 
the original singularty distribution. 

The boundary condition on each body gives no flow 
across the boundary which describes the rigid body. In our 
case, this is a body of revolution. At each location along 
the body there is a radial distance at which the velocity 
(radial) must be zero. If the cross flows were uniform, 
the motion of an infinite cylinder could be represented 
exactly in the fluid by placing the appropriate doublet 
distribution on the axis, The flow patter around a slender 
body of revolution can be approximated locally by an infi- 
nite cylinder because the change in diameter would be small 
for a change in the axial direction. For the distances 
between bodies of interest, it was estimated that the uni- 
form flow assumption for any one segment of the body would 
be nearly correct, The dipoles were sized with flow of the 
value experienced on the body axis as a measure of the 
uniform flow strength. For the actual flow conditions, the 
boundary location is not calculated. It is assumed that 
the deviation from the exact position would be small. The 
dipole distribution is sized by the following rule, See 


Appendix C for derivation, 


SD, (x = -2 * V(x) * AREA (x) / av 
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WHERE AREA = CROSS SECTIONAL AREA 
V 


tt) 


CROSS Flow VELOCTTY IN FT/sec, 


If the perturbed flow is known, the changes in the 
Singularity distribution can be readily determined, For 
the case of two bodies near each other, the change in the 
Singularity distribution in one of the bodies is reflected 
in a slight change in the perturbed flow field experienced 
by the others. What is desired is a simultaneous solution 
to the problem. This is achieved by performing several 
meeratvons of thesatoaw, Taking farst one, than the other 
body, change its singularity distribution to counteract 
the perturbed flow, and then determine the resulting 
affect on the flow field at the other body. Changes in 
the perturbed flow were insignificant after about three to 
five iterations except when the bodies were etremely close 
and then it took only about ten iterations to achieve con- 


vergence to the solution. 


Pen intereculenm Force. Calculations 


Mieminterdculom forces actins On a rigid body moving 
in an ideal fluid without viscosity can be determined if 


the fluid velocity field can be represented by a potential 
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function, In principle, the method for determining a force 
on any body would involve the calculation of the pressure 

on the body surface followed by the integration to obtain 
the total force. This was first accomplished by either 
Munk or Lagally for bodies generated with singularities. 
Lagally noted that for a bocy placed in a steady stream 

this integration took a form that would be analogous to a 
derivation which treats a force as acting between the sin- 
gularities. The apparent force between a pair of Singulari- 
ties is determined in three dimensions by the following 


formula, 
Fe ene (eqaenen / 2 


Where a poSitive force indicates 
an attraction, 


The total force would be determined by adding all 
pair contributions. The force acting on a body which is 
represented by a continuous source distribution along its 


axis has the following integral form. 


bb mx mi(y') 
a dxdx’ 
Pei Lf priry 


Where [a,b] and f[a',b'] indicate 
the line segments with singularities. 
The results of Lagally were expanded by Cummins to 
include the forces between doublet singularities for the 


general time dependent case. His derivation involved the 
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following assumptions, 


ie unesvelocity Tield is irrotational 
and has a velocity potential % (x,y,z,t) 


2. if the body were not present, the 
stream would have a velocity potential ¢, 
which we call the potential of the "undis- 
turbed stream", 

3. There are no Singularities of the 
undisturbed stream in the region occupied 
by the body. 

4, The boundary condition at the sur- 
face of the body is satisfied by super- 
imposing a system of singularities upon 
the undisturbed stream, such singularities 
falling within the region which the body 


would occupy. The potential of the system 
of singularities is designated by gp. Then 


$= 6 + Br 


Cummins found that the forces and moments calculated 
were of the same form as expressed in the Lagally deriva- 
tion, but with additional terms, He expressed the Lagally 
force term as a product between local Singularities and 
the total flow perturbation rather than as the sum of in- 
dividual contributions from each singularity, The two 
representations for the force are identical, Zucker gives 
a summary of the results and indicates the contribution 
by Landweber and Yih in substantiating Cummins' results 
and evaluating the moment term resulting from the change 
in the fluid flow with time. 

The force on the body will be as follows for Cummins' 


assumptions, 


> 





ae e¥ a Coke Sere bs 


e Li, ah] +m, +M.+ Ms 


Sh 


WHERE 
VY = Volume of the body 


> Is the translational velocity of the body 


= Is the location of the centroid of the 
body relative to the origin of the moving 
system, 


F = The force if the body were not rotating 
and the undisturbed stream were steady 
(Lagally Force) 


F,7Force due to change of the flow with time 


F,-Force arising from body rotations 


M,=Moment if the dlow were steady (Lagally) 
M, = Moment due to the flow changing with 
time 


M,* Moment effect due to rotation of the 
body 


FOR SITUATIONS DESCRIBED BY ONLY SOURCES AND DOUBLETS, 


Ps AND NM, ARE ZERO. 
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FOR A SOURCE OF STRENGTH 
ee 
ee een ea 


defines the fluid velocity in the 
undisturbed stream 
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WHERE Pz JWWAAD DIRECTED NORMAL VECTOR 


This term was oval Nata Landwebber and Yih, 





The force calculations programmed on the computer 
used at least one of the forms discussed previously, As 
more complicated and detailed models were developed the 
form of the force calculation was selected to improve com- 
putational efficiency. For the Rankine Ovoid and ellipsoid 
which do not contain cross flow corrections, the simplest 
statement of Lagally's force was used. When it became 
necessary to introduce dipoles to correct for the cross 
flow, the final form of Cummins was used, The final 
slender body model takes into account the terms that are 
represented by F, and M, . This represents the non-zero 
terms of the interaction force for steady parallel motion 
of the bodies, 

The above equations relate to a three dimensional 
flow problem. The two dimensional case is Similar with 
the Lagally force being of the following form: 

Fi, = aire m2 fp 
Havelock developed in Reference 7 the force and moments 
acting on a two dimensional body by using contour integra- 


Gon. 


mone sOURCE DISTRIBUTION 
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WHERE 


POV t ie ow ebrerers to the given distribution 
in the liquid, and "r" to the image system 
within the surface of the body. 

2, 


Summation extends over the external and inter- 
nal sources taken in pairs, 


Eon DIPOLE DISTRIBUTION 


| 


a y = es iM M. cos(%,.+Hs-3O,, 
an Gy TS elie) 


~< 
0 
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7 
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MM, sv Gr tas ~30,, 
447reZ Me Ms )/p3 


(r,5) (3) 


Where; — the doublet is of moment M and makes 


an angle X with the x axis (ship). 


And: ©, being the angle between x’ and the vector 


drawn from the doublet (r) to the doublet (s). 


The total moments results from moments of forces already 


determined plus an additional term, 


Mi = air C 2.M.Ms Siv (x, 4X; - 28,, ) Jae 
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2.3 Computation 


The earlier discussion of the slender body model gave 
the general steps involved for determining the interaction 
forces, In this procedure there were two major areas 
requiring considerable effort, First, it was necessary that 
derivatives of the flow potentials be derived for the force 
calculations. These potentials and their derivatives were 
defined in terms of a moving axis system of the ship in 
which the forces were acting. The necessary transformations 
and derivations are developed in Appendix B, 

The second area involved the determination of the total 
force acting on the ship. The slender bocy was defined by 
a continuous distribution of singularities, The resulting 
flows in the fluid were determined by integrating the con- 
tributions from the singularity distribution, Also the 
forces were determined by integrating the individual force 
contributions. In both cases, the numerical integration 
used values of the integrand at the same location on the 
body axis. These coincided with the stations at which 
the sectional area curve had been defined, The integration 
was performed in each case using Simpson's Rule with an 
odd number of stations (21), evenly spaced along the body 


axis, 


he 





CHAIvieie 101. 
COMPARISON OF RESULTS 


371) eayvlor Model. Tests 


The Rankine Ovoid and ellipsoid models were developed 
and checked with the available results then used to model 
the forces and moments reported by Taylor, Figure 7 de- 
picts the results. These predictions cdiffer in two ways 
with Taylor, First, they do not predict the secondary re- 
versal in the moment that is evident at longitudinal se- 
parations of .6*L. Secondly, the largest moment is pre- 
dicted for longitudinal separations above 1.0.x L. This 
variation can be attributed to the deformation of the boun- 
daries of the bodies as they approach each othr. If the 
boundary was maintained in the original position, the flow 
would have been restricted, thus increasing the fluid 
velocity and contributing to a lower pressure on the bodies. 
This reduction in pressure would contribute a moment that 
would cause a reversal in sign. The inability of the pre- 
vious theoretical models to explain the reversal in sign 
of the moment curve, led to the development of a model which 
incorporated cross flow corrections. The correction for 
cross flow indicates the reversal in the moment can be 
modeled. However, the magnitudes were not in full agree- 


ment with those obtained by Taylor. Computations were made 
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at higher speeds with no variation in the results, This 
was expected, the plot involved the non-dimensional moment 
po ometsient. For closer distances, the predicted ampli- 
tude of the moments increased only slightly. (Yhe model 
tests were conducted at a very small separation. ) 

The dimensions of the models used in Taylor's experi- 
ments were inputs to the slender body model. Calculations 
were made for the same cases as described by Taylor, The 
following curves show the results plotted for comparison 
purposes, 

Additional curves are drawn which show the forces pre- 
dicted for the Taylor full bodied models, as the lateral 
separation is increased, The reversal in moment diminishes 
as the ships separate because there is no longer a restric- 
tion in the flow to cause a pressure drop between the ships. 
Beyond a distance of 2 beams, the results agreed with those 
obtained with a model which did not incorporate cross flow 
corrections. Only beyond a separation of 4 beam widths 
interaction effects were negligable. 

The model tests conducted by Taylor, and the slender 
body predictions differ in several ways. Some of the more 
Significant assumptions used for the slender body were: 

1, The potential field was assumed most 
important (No free surface) 


2. The ships were represented by a body 
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of revolution, (Lessen effective beam) 
3. Squat and sinkage were not incorpo- 

rated, 
Other differences in the results may be attributed to either 
asemmetry in the models or inaccuracies in the measurement 
techniques, The apparatus used by Taylor was not indicated, 
When the tests were conducted, the resistance measurements 
used a system of weights, pulleys and levers, Taylor even 


commented on the problem of conducting the model tests, 
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3.2 Comparison of Results with Newman 


In order that a comparison be made between the model 
developed and the results of Newman. The sample problem 
in Newman's paper was used, The sectional area curve for 
the ships moving in parallel were determined using the 
radius formula presented earlier. Two ships moving in 
parallel were defined with the same geometry. The initial 
calculation assumed a length of 600 feet, a beam of 100 
feet and a speed of 10 knots. The lateral displacement 
was determined and plotted according to the non-dimensiona- 
lization used by Newman, 

Note that Z,./ro is equivalent to the separation dis- 
tance between the axis of the two bodies divided by their 
beam dimension, Other beam to to length ratios were cal- 
culated by varying the length of the bodies while main- 
taining the beams constant. These included the L/B of 6, 
Mee ormoo, The tetal force on the body of revolution 
determined using the model. was non-dimensionalized by the 
maximum value determined by Newman. The following graphs 
depict the results developed in this thesis, 

The first graph shows the results for a L/B of 6 for 
the various theoretical models, Calculations were then 
made for the resulting forces with and without the cross- 
flow correction which used the dipole distribution. It 


was expected that the results of Newman and the proposed 
1 
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model would approach each other in the limit of high L/B 
ratios and for very slender bodies because the singularity 
distributions would be nearly identical and lie very close 
to that used by Newman. 

The next graph shows the change in the non-dimensiona- 
lized maximum force as a function of the length to beam 
ratio, Note that the singularity model approaches the New- 
man result in the limit. On the graph, the error of the 
computation is indicated. The solution with the dipole 
correction seems to converge initially and then for very 
high ratios diverge. This tends to indicate that the di- 
pole correction may be useful for all ratios of interest, 

The error introduced during integration was easy to 
approximate for the axial distribution of singularities 
because the fourth derivative of the integrand should be 
exactly zero. However, for the case with the dipole cross 
flow correction this may not be the case and it is expected 
that the introduced error would be larger. This may ex- 
plain the divergence of the solution for high L/B (large 
lengths). See Appendix E. for derivation of error estimate. 

Another graph shows more directly the variation in 
the results for the various length to beam ratios, In this 
graph the shaded area indicates the effect of the dipole 
correction, 

The last graph is a compilation from various sources 


of interaction forces. Tae plotted in the same 
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non-dimensional terms. For the ships, the forces were 
doubled to approximate the forces that would be observed 
between bodies of revolution. The non-dimensionalization 
with respect to the maximum force predicted by Newman is 
consistant with the other graphs, FPlotted also are the 
solutions predicted for a body of revolution with L/8 
ratios of 6 and 10 with the geometry used by Newman. The 
plot indicates that the model tests agree more closely to 
the results predicted by the proposed model with the results 
of the tests by Newton showing the greatest variation, If 
the results had been mislabeled the results would have been 
more consistant with the proposed model predictions. 

The basic differences in the proposed model and that 
of Newman are, first, the surface condition on the body is 
accounted for by displacing the singularities from the 
centerline in the case of Newman, where the proposed model 
uses a dipole distribution, It is obvious that the pro- 
posed model will satisfy the surface condition better when 
the induced velocity is fairly uniform over the body pro- 
file, From the data presented, the relatively small L/B 
and the separation distances of interest justify the use 
of dipole corrections. A second difference in the models 
is in the evaluation of the integral solutions, Newman 
uses an approximation to the solution which emphasizes 
the need to restrict the solution to the slender bodies, 


whereas this model uses pee ere tion. 
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CONCLUSIONS _ 


The major effort of the thesis has been to provide a 
theoretical model for the prediction of forces between 
ships. A computer program has been developed which carries 
out the computations for ships in rectilinear motion using 
the potential flow model described. The ellipsoid model 
proposed by Havelock was unable to predict the form of the 
results as indicated by empirical model tests. These 
tests demonstrated a reversal in the moment force as the 
ship separation decreased. This led to the development of 
the slender body model with cross flow corrections. The 
programs which were developed for ships in rectilinear mo- 
tion ignore time dependent forces. The data presented 
here represents the interactions of ships in deep water, 
on parallel courses and with each ship at the same speed, 
Orientation, other than parallel, would have resulted in 
a change in the relative positions of the ships with time 
and would involve other force terms as derived by Cummins, 
but were not programmed here. 

The investigation of past efforts indicate a lack of 
good model or full scale data to test a theoretical result. 


The differences experienced in the analysis of Taylor's 
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results could not be fully explained without more empirical 
data. This can only be obtained with ship models because 
of the parallel course requirements, 

In meeting the objectives of the thesis, the axial 
distribution of singularities and the force calculation 
techniques provide additional inputs to a general trajec- 
tory model, When the ship approaches other bodies or 
boundaries in the fluid, this work will help to predict 
the interactive forces so that more inclusive trajectory 
calculations can be made, This is applicable to the deter- 
mination of the factors involved in ship collisions, the 
prediction of the control parameters for at sea replenish- 
ments, or for future ship design characteristics, 

The performance of model tests may indicate results 
that differ from Taylor's. The design of model experi- 
ments must take into account the large forces introduced 
by any misalignment of the models, Also, the wall effects 
may introduce erroneous readings, The computer program 
has been written to include more than two ships in proxi- 
mity. It could be used to estimate the wall effects in 
the tow tank experiments by introducing the first few image 
distributions to approximate the wall boundary conditions, 
The following figure shows the orientation, It will be 
important that some tests be repeated with the models in 
a reverse geometry to insure accurate readings. This would 
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FIGURE 14 


THEORETICAL CALCULATION OF WALL EFFECTS 
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give a measure of the wall effects and asymmetries in the 
models, 

The program can be extended to cover shallow water 
effects, For very shallow water, a two dimensional model 
could be developed by only a few changes in the present 
program, These would involve the change in sizing the 
Singularity distributions and in the Lagally force calcu- 
lation, For other shallow water conditions, the boundary 
conditions could be satisfied by a series of image distri- 


butions as shown in Figure /5. 
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FIGURE 15 
SHALLOW WATER APPROXIMATION 
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CHAPTER. 


RECOMMENDATIONS 


The work accomplished in the thesis points out the 
need for continued study of the interaction between ships. 
The following list indicates possible areas for future in- 


vestigation. 


1. Perform model experiments to obtain data to 
verify Taylor's results, including shallow water effects 
as described, 

2. Include remaining force terms.in order to calculate 
for unsteady ship motions. 

3. Incorporate program in trajectory model. 

4, Determine stable positions alongside, and test 
with operating experience. 

5. Simulate refueling maneuvers and compare with 
operator experience. 

6. As a final test, attempt to reconstruct published 
collisions involving ship interactions, 

7. Investigate feasibility of extending the theoreti- 
cal model to include free surface condition by distributing 
Singularities as used by Havelock and developing the re- 
quired theorems to evaluate the interaction forces as accom- 


plished by R. D. Zucker for vortex sheets, 


te 





O. 


10, 


BIBLIOGRAPHY 


G. K. Batchelor, An Introduction to Fluid Dynamics, 
Cambridge University Press, Cambridge, England, 1967. 


C. N, Calvano, "An Investigation of the Stability of a 
System of Two Ships Employing Automatic Control While 
on ParalleJ. Courses in Close Proximity", M.1.T. Thesis, 


omen 


W. E. Cummins, "The Forces and Moments Acting ona 
Body Moving in an Arbitrary Potential Stream", 
David Taylor Model Basin. Report No. 780, Washing- 
Pomme OC. , 1953. 


R,. S, Garthune, et. al., "The Ferformance of Model 
Ships in Restricted Channels in Relation to the 
Design of a Ship Canal", David Taylor Model Basin 
Report No. 601, Washington, D. C., 1948, 


R, W. L. Gawn, "The Admirality Experiment Works, 
Haslar"”, Quarterly Transactions of the Institute 

of Naval Architects, Vol. 97, No. 1, London, England, 
January, 1955. 


M. Gertler, "A Reanalysis of the Original Test Data 
For the Taylor Standard Series", David Taylor Model 
Basin Report No. 806, Washington, D. C., March, 


1954, 


T. H. Havelock, "The Vertical Force on a Cylinder 
Submerged in a Uniform Stream", The Collected Papers 
of Sir Thomas Havelock on Hydrodynamics, Edited by 


C. Wigley, Office of Naval Research, Department of 
the Navy, ONR/ACR-103, 1963, pp. 297-303. 


, Discussion to "Interaction Between Ships”, 
by A. M, Robb, Transactions of the Institute of Naval 
Architects, Vol. 91, London, England, 1949, pp. 336- 
Be. | 


Sir H. Lamb, Hydrodynamics, 6th Edition, Dover Publi- 
cations, New York, New York, 1932. 

H. T. Larson and L. M. Tichvinsky, "Replenishment at 
sea; Experiments and Comments", University of Cali- 
fornia Institute of Engineering Research Report No. 
154-1, Berkeley, California, June 1, 1960. 


58 





ee 


eae 


13. 


14, 


15. 


HO, 


Lie 


J. N. Newman, "The Force and Moment on a Slender 
Body of Revolution Moving Near a Wall", David Taylor 
Model Basin Report No. 2127, Washington, D. C., 1965. 


R,. N. Newton, "Some Notes on Interaction Effects Be- 
tween Ships Close Aboard in Deep Water", First 
Symposium on Ship Maneuverability, David Taylor 
Model Basin Report No. 1461, Washington, D. C., 
October, 1960, pp. 1-24, 


L, Prandtl and 0. G. Tietjens, Fundamentals of Hydro 
and Aero Mechanics, McGraw Hill Book Company, Inc., 
New York, New York, 1934. 





B, iL. Silverstein, "Linearized Theory of the Inter- 
action of Ships", University of California Institute 
of Engineering Research Report 82-3, Berkeley, Cali- 
homta, Mave 15... 1957. 


V. L. Streeter, Fluid Dynamics, McGraw-Hill Book 
Company, Inc., New York, New York, 1948, 


D. W. Taylor, "Some Model Experiments on Suction of 
Vessels", Transaction of the Society of Naval Archi- 
tects and Marine Engineers, Vol. 17, 1909, pp. 2-22. 


R. D. Zucker, “Lagally's Theorem and the Lifting 


Body Problem", Journal of Ship Research, Vol. 14, 
June, 1970, pp. 135-141. 


De 








APPENDICES 


60 








SsaAc (I,J) 


SAC 


DSAC (I,J) 
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LENGT (TI) 
BEAM (TI) 
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CM (TI) 
Drom (1) 
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SM (I,J) 


Sp (1,J) 


NN (TI) 


KK (I) 
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Gin. ) 


(I,J) 


APPENDIX A 
Nondimensional sectional area at Station 
Otis orth oaeel 
Sectional area at Station J of Ship I 


Derivative of sectional aréa at Station J 
Gi ship 1 


Displacement of Ship I (Not used) 
Length of Ship I 

Beam of Ship I 

Velocity (retilinear) of Ship I 


Prismatic coefficient of Ship I 
Midship coefficent of Ship I 


Distance from midship where R. 0. singu- 
larities are placed to generate Ship I 


Initial (no interaction) source strength 
at position J on the axis of Ship I 


source strength at position J on the axis 
of Ship I 


Dipole strength at position J on the axis 
of olip es! 


Number of stations used to input non- 
dimensional sectional area curve and to 
perform calculations 


Index for Ship I to indicate model used 
to calculate forces (not used) 


Mapes tiOn.of sonap J in. the fluid 
Mmpoatcloneot ofp 1 in the fluid 


Heading angle of Ship I in the fluid 
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PI = The constant 


RON = Fluid density 

M = Index. (not used) 

MM = Index (not used) 

N = Number of ships interacting 
(Program can handle up to 5 with space allo- 

cation) 

Meas) = Compencme of force ian the X direction at 
station J 

Pit) = Component ot force in the Y direction at 
Station J 

RF1(J) = Moment at Station J on Ship I 

XF2 = Component of force in the X direction of Ship 2 

J = Index for indicating position along Ship I 

L = Index for indicating position along Ship K 

X1i(J) = Longitudinal distance to Station J on Ship I 

Vick = Lonamgndinal distance to Station & on Ship K 

I = Index to indicate first ship of the pair for 
force calculation 

K = Index to indicate second ship of the pair for 
force calculation 

AFX(I) = Total force on Ship (I) in X direction (Out- 
put on non-dimensionalized ED-form) 

AG) Tota la toree on Ship (1) in Y direction 

AM(I) = Total moment on Ship I 

ACM(I) = Total non-dimensionalized moment for Ship I 
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DGNZ (J) 
QRZ (J) 
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DXI (J) 


DXK (L) 


AEX (T) 


It 


(jtclenlov az POSition J of ohip 2 (not 
used ) 


Total axial flow at Position J of Ship I 


Total transverse flow at Position J of 
Ship 1 

Total axial contribution of flow gradient 
at location J of Ship 1 

Total transverse contribution of flow 
gradient at Location J of Ship i 


Axial flow contribution at J on Ship I 
due to singularity at L on Ship K 


Transverse flow contribution 
Transverse gradient contribution 
Flow contribution 

Axial gradient contribution 


Station Spacing on Ship I between Station 
J and Jti 


station spacing on Ship K between L and 
L+1 


Focal distance of ellipse measured from 
midship position 
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APPENDIX B 


Derivation of Axis System and 


Transformations and Potential derivatives 
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MOVING AXIS SYSTEM 
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The distributions used to model the Taylor full ship 
results are as follows, The situation when the models 
are abeam with a 3.9 ft. separation is represented (end 
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APPENDIX F 


Presentation of Ellipsoid and 
Rankine Ovoid Results (Compare with original work) 
The results are the same except for 
a scale change. Havelock's scale should be divi- 
ded by 10, where Silverstien's should be divided 
vy (47r)°, or 158, 
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APPENDIX G 


COMPUTATION E°FICIENCY 


The programs written without cross flow corrections 
take only a fraction of a second to evaluate interaction 
forces for each individual ship position. The slender 
body model with cross flow corrections takes various times 
for execution depending on the number of iterations of 
the flows required. For the samples run, two ships with 
twenty one stations each, took three seconds for five 
iterations increased to ten, the time doubled to six sec- 


onds for each ship position calculated. 
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APPENDIX H 


Applications Program Series 
AP-21 revision 2 


NAME : 


DESCRIPTION : 


METHOD: 


USAGE: 


MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
INFORMATION PROCESSING CENTER 


September 28, 1973 
ZEROIN, a double precision FORTRAN IV subroutine 
to solve a system of simultaneous, nonlinear equa- 


tions 


ZEROIN computes a solution vector X of a set of N 
simultaneous, nonlinear equations: 


F,(x)=0, 
Pa (X)=0. 
X is obtained by an iterative method: 

ey eS x4 pik) 
where the vector p!k) = gk Yap x) y | with 
ac a representing the vector of function 
values at the point te and yon representing 


an approximation to the Jacobian at x () The 
following test is used to determine ees cen 


litem) _ x 8) | 2 < eps«||x(K*2) | )2 


where EPS is a user-supplied tolerance. All real. 
arithmetic in ZEROIN is double-precision. 


The call for ZEROIN is: 
CALL ZEROIN(N,X,R,EPS,ICV,EVAL,IRITE) 


Input values are: 


N the order of the system of Bail 
2 < N Ss. ZU 

Xx a one-dimensional array initially con- 
taining an estimate of the solution 
vector. 
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Pee 


EPS 


RITE 


EVAL 


Output values 
X 


R 


EPS 


ICV 


Page 2 
an estimate for ||solution vector-x| | 


a tolerance to be used in the test 
for convergence, usually 90.5E-12. 


a specification for intermediate 
printing. If IRITE = 9, ZEROIN does 
NOspranting. li VIRITE = 1,..-2ZEPOIN 
prints, at each iteration, the current 
value of R, the current approximation 
to the solution, the number of times 
EVAL (below) has been called, and the 
current function values F(X). ° The 
printed output from the final itera- 
tion contains, in addition, the output 
value of ICV. 


is the name of a user-supplied 
SUBROUTINE subprogram that evaluates 
the functions F at each iteration. 

The call for EVAL is CALL EVAL(X,F) 
where X is the input to EVAL, and F 
contains the function values corres- 
ponding to X. xX and F must be singly- 
dimensioned, double-precision vectors. 


will be: 


the solution vector. 
k+l k k+l 
sees | 120i jx aie 


where xe) is the solution vector. 

When ZEROIN terminates normally, the out- 
put value of R will be less than the in- 
itial value of EPS. In‘case of abnormal 
termination, R can be used to determine 
what degree of convergence was obtained 
in the output value of X. 


| |rcx) |]? where F(X) is the vector of 
function values corresponding to the 
solution vector X. 


one of the following indications: 


icv = 0 if convergence was obtained, but 
|JF¢(x) | ]* > 0.5E-10 


Icv = 1 if convergence was obtained, and 
|JE(x) ||? < 0.5E-10 


Icv = 2 if convergence was not obtained 
within M iterations, where 
91 M = MIN(10*N,50). 





Rea el 2 


CODING 
CONSIDERATIONS: 


AVAILABILITY: 


CORE 
REQUIREMENT: 


REFERENCES: 


AUTHOR: 


Appropriate declarations are: 


REAL* 8 X(N), R, EPS, EVAL 
INTEGER* 4 N, ICV, IRITE 
EXTERNAL EVAL 


1. The input value of R may be chosen to be 
||X||/10.0, where X is the input vector of es- 
timates, ig the solutions are all about the 
same magnitude. 


2. The output value of R can be used to determine 
degree of convergence in the case ICV=2. 


3. ZEROIN is more likely to converge if the equa- 
tions are scaled so all solutions are of the 
Same magnitude, and if the initial estimates 
are good. Poor initial estimates may produce 
divide check interrupts during execution of 
ZEROIN. For details on the choice of initial 
estimates, refer to the first reference below. 


ZEROIN is available in load module form in the li- 
brary SYSS.MATHLIB.SUBR and in source form in 
SYSS.MATHLIB.SOURCE. For further information about 
the use of these libraries, refer to: AP-69 or con- 
sult the Programming Assistants in Room 39-219. 


On the 370/165 ZEROIN uses 9,730 bytes of core stor- 
age (hexadecimal 2602). 


Mancino, O.G. “Resolution by Iteration of Some Non- 


ee 


Computing Machinery, Vol. 14, No. 2 (Npril, 1567), 
pages 341-350. 


Robinson, S.M. "Interpolative Solution of Systems 
of Nonlinear Equations," SIAM Journal on Numerical 
Analysis, Vol. 3, No. 4 (December, 1966), “page 650, 


ZEROIN was written by G.W. Westley, Computing Tech- 
nology Center, Union Carbide Corporation, Nuclear 
Division, Oak Ridge, Tennessee. 
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Computer Flow Chart 
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